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particles are tightly coupled to the flow, then the collecting probe receives 
information in electrical form about the fluid velocity. Besides being a 
direct velocity measurement, the probe can be of streamlined shape and very 
rugged, and it is insensitive to temperature differences in the flow. The 
charged particles can be chosen so as to give the desired range of frequency 
response. The instrument is sensitive to impurities in the flow but in 
principle the effect may be calibrated out. 

An EGD spectral anemometer has been tested using water droplets charged 
as they condense through a corona discharge. The anemometer consists of a 
short-circuited EGD generator; this is the first time that such an application 
has been reported and a Navy patent has been awarded for this device. The 
feasibility of the anemometer depends on the magnitude of the charged particle 
slip. The coupling between the flow field and the electric field depends on 
the frequency range of the velocity fluctuations as well as on the charge, 
size, and mass of the particles. The work treats the mass as a first order 
perturbation, in the particle momentum equation. The unsteady flow studied is 
a region in the wake of a cylinder where charged water droplets are injected 
into an air flow. From this work the size of the charged particles is calcu- 
lated and a correspondence is found between the mean particle size injected 
and the observed frequency response of the anemometer tested. The EGD spectral 
anemometer data are compared to hot-wire data taken at the same location; some 
results compare within + 2 dbs up to a frequency of 630 Hz whereas others are 
within + 6 dbs up to 20 KHz. The parameters which varied are the steam 
pressure and temperature, and the corona characteristics. Predictions of 
frequency response as a function of charged particle size are given. 
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I. INTRODUCTION 



Electrofluiddynamic devices operate on the principle that a 
nonconducting fluid can interact with an electric field when unipolar 
charges are introduced into the fluid. If this interaction does not 
disturb the fluid appreciably and if the charged particles are linked 
intimately with the flow, then the arrival of these particles at a 
suitable collector brings information about characteristics of the flow. 
This work deals exclusively with gases and the term electrogasdynamic 
(EGD) anemometer will be used to denote the device under consideration. 

The feasibility of the anemometer hinges on the proper choice of 
charged particles. These must couple the gas to the electric field 
through individual collisions if the particles are small enough or 
through fluid dynamic drag and pressure forces, if they are large enough. 
Molecular ions are too mobile to be useful^"; therefore, charged particles 
of size much larger than the gas molecules must be introduced. The 
ultimate size of these particles is restricted by their ability to 
follow the turbulent fluctuations and by the fact that, depending on the 
method of charging, a large particle may have a sizeable mobility (as 
will be seen in Section III). In this paper, the effect that particle 
size and charge have on the coupling between the fluid and the electric 
field will be discussed. 

It may be mentioned in passing that the definition of particle slip 
in unsteady flow is not only crucial for the anemometer but also for the 
EGD-power generator. Since turbulence is a means of enhancing the break - 

2 3 1 4 

down strength in a gas ’ and since generators ’ rely on the mixing of 
fluids both to increase the pressure drop in the channel and to insulate, 
there must exist regions of unsteadiness in -the flow. Hence, a small 
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slip parameter for steady flow may misrepresent the ability of a charged 
particle to interact with the flow. 

II. EC-D - Anemometer 

When a charged particle is tightly coupled to the flow, it is 
essentially a flow "tracer." To pick up such a tracer one simply places 
a metallic probe which provides the neutralizing current. The injected 
charges need a certain time to equilibrate with the flow but there- 
after they drain no energy from the flow. The particles 

are neutralized as they arrive with the surrounding gas and are not 
forced to the probe by an external, field hence, in principle, the probe does not 
distrub the nature of the flow. 

Because the collection area is not necessarily the probe area'' - and be- 
cause of an uncertainty in the size and charge of the particles, 
it is expedient to measure the spectral characteristics of the flow 
independently with a hot-wire anemometer and then to compare the results 
with those of the EGD -anemometer . This permits the optimization of the 
operating parameters of the EGD -anemometer and allows the inference of 
charge and size distributions for the charged particles. 

The hot-wire anemometer, although a heat transfer measurement, has 
emerged as perhaps the most useful device for measurements in unsteady 
flows. Needless to say, advances in fluid mechanics hinge on reliable 
measurements, and hot-wire data are an accepted standard for these 
measurements. The EGD device can complement the hot-wire anemometer. 

An EGD -anemometer consists of the following: an injector or source 

for the charged particles, a conversion region where the particles 
equilibrate with the flow, and a collector for the charges. The collector 
unit corresponds to the anemometer probe. The EGD-probe is very rugged, 
allowing the probing of intense turbulent fields. It is simple and 
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inexpensive, and therefore can be easily manufactured. Since very small 
charged particles can react to high frequency oscillations, the frequency 
response of the EGD- probe can be as high as the best hot-wire (at no 
compromise in mechanical strength). In addition, signals from the EGD- 
probe are electrical and need little or no further amplification. Effects 
corresponding to the thermal inertia of hot-wires are not present. There 
is no need for a linearizer since the EGD measurement is already linear, i.e 
a direct measurement. Absent are the complications that free convection 
currents bring to the hot-wire zero setting. In the case of two-phase flows 
the EGD-probe needs no separate calibration and will withstand the impact 
of charge droplets or solid particles whereas the fragile hot-wire will 
readily break under these conditions. In the case of two mixing jets, 
there will invariably be a temperature difference and this, however slight, 
will affect the hot-wire since it is sensitive to temperature differences. 
The geometry of the EGD probe is more adaptable than that of the hot-wire 
to unusual probing situations. 

The main disadvantage of the EGD scheme for anemometry is the require- 
ment for a source of particles of optimum size and charge. If the particles 
are created through a corona discharge then there is a need for a high 
voltage supply (with the associated safety precautions). Also, the 
injector/ collector unit must have a separation of the order of millimeters 
and must be traversable all of which pose several design problems. However, 
all these problems can be overcome and the use of an EGD-generator scheme 
as a spectral anemometer looks quite feasible. 

It should be pointed out that the use of charges in a flow for 

anemometry is not entirely a new idea. Hinzd^ describes electric discharge 

anemometers which consist typically of a glow or corona discharge in the 

£ 

flowing medium. More recently, Nygaard has worked on a type of corona- 

discharge anemometer. A mass flow meter device based on an ion- displacement 

7 

concept of Professor Durbin is being marketed by Thermo- Systems, Inc. 

4 



None of the reported schemes, however, have concentrated on the spectral 
response of the instrument and, in particular, in matching the charged 
particles to the desired value of the slip at a given frequency range. 
III. Concepts from Unsteady Flow 

The ability of a charged particle to couple its surrounding gas to 
the electric field is described by its slip. The drift velocity with 
which a charged particle moves in an electric field is given by the 
product of its mobility (p,) times the electric field (E) 

v D = p.E (1) 

and the slip (S) for steady flow is defined as 




V = fluid velocity 

q 

Here one wishes to make S as small as possible. 

So, knowing the electric field and the fluid velocity, the slip can 
be calculated provided that the mobility is also know. Unfortunately, 
deducing the mobility is not straightforward since it depends on the size 
of the particle and on how it was charged'*'. If one plots the drift 
velocity defined in Equation 1 for the maximum field (breakdown) as a 
function of the radius for spherical water droplets in air at standard 
conditions, one gets the curves shown in Figure 1. Singly charged 
particles have a decreasing (iE^ -velocity with radius because the drag is 
increasing while the electrical force remains constant. The case for 
particles charged to the Rayleight limit represents an upper bound in the 
charging scheme because more charges exceed the stability of the droplet 
at any given radius. Corona field charging is an intermediate case 
between the other two modes of charging. • Since the speed of sound is 
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FIGURE |. MAXIMUM DRIFT DUE TO ELECTRIC FIELD FOR 
WATER DROPLETS IN AIR AT I -ATM AND 20° C 
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representative of flow velocities of practical interest. Figure 1 gives 
information about the range of radii and the charging scheme that would 
be required for a given value of the slip. 

When the flow is unsteady, the question arises as to what happens 
to the conventional slip parameter. Clearly, inertial effects of the 
particles will come into play. In Appendix A, the momentum equation for 
a particle is treated. There a one -dimensional, spectral analysis is 
given which shows that the concept of the slip -parameter may be extended 
to unsteady flow. Written in a spectral form, the maximum value of the 
slip becomes 

p,K 

S(co) I = TT7 v + CUT (2) 

v ' 'max V(u>) v 

where 

cu = frequency 

V(cu) = RMS spectral intensity 

t = e-folding time for the particle velocity to approach 
steady state value 

(the equilibration time) is a function of the particle radius and 
of the viscosity of the fluid. If turbulence or any unsteady flow is to 
have an effect on the charges, then the slip parameter must be small for 
all the dominant frequencies (cu^) of the flow, i.e., 

(J.E 

+ \ T v << 1,0 

Actually, it is only necessary to look at the highest frequency of 

interest since for all smaller values of cu, the inequality above will be 

satisfied provided V(cu^) retains nearly the same magnitude. 

Figure 2 shows two curves bracketing a region in which the 

-3 -2 

contribution of cut_^ to Equation 2 is between 10 and 10 . This range 

represents a criterion of small slip. The contribution of the first 
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FIGURE 2. OPTIMUM FREQUENCY RESPONSE 

AS A FUNCTION OF THE RADIUS 
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term in the right-hand side of Equation 2 may or may not be as low as 



with the charge of the particles (see Figure l). 

In Figures 1 and 2, water droplets have been described. 

For less dense particles or for hollow, spherical solids the curves move 

to the right (i.e., towards larger radii) which is a desirable effect 

from practical considerations. It is conceivable that one size of 

particle will be optimum for many applications. 

The measurement of particle size, or more accurately stated of 

particle size distribution, is very difficult in the range of interest 

9 10 

of EGD-anemometry 7 ’ . The water droplets can only be measured in situ 

and direct techniques such as laser holography'*''*' can presently not 

-6 

measure below 50 x 10 m*. In the course of this work an indirect 

method of measurement was developed which consists of comparing 
the spectrum of the current collected to that of the hot-wire 
spectrum. For the proper particle size and charge distribution, the 
agreement between the two measurements is high when the flow has 
not been significantly disturbed by the particles. 

The current collected may be written as 



_2 

10 but one would endeavor to pick a small enough value by manipulating 
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(3) 



where 



R = radius of the particle. 




00 




and 



q(R) = charge of the particle 
Vp(u},R) = velocity of the particle 
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Dividing Equation 3 by v(oj), the flow velocity, which is independent 
of the particle radius, one gets 



I (go) 

g 

v(tu) 



-Jo 



,Kr) 



q(R)n(R)dR 



(4) 



It is shown in Appendix B that the ratio v (go ,R)/v(oo) may be approximated 
by a Heaviside -type of function, 



VpKR) 

v(«>) 



= 1 
= 0 



T „ * V‘ 



O) 



T > l/lO 
V ' 

This defines R’(u)), the radius at which equals w ^ and Equation 4 
becomes 

I (<u) R»(go) 



g 

v(«>) 



J o q(R)n(R)dR 



where for spherical particles 

. . 

= viscosity of fluid 
p = particle density 

Since one can obtain I (oi) and v(cu) experimentally, and since one 

§ 

can surmise what q(R) must be, one can deduce the distribtuion n(R) by 
solving the implicit equation above by successive approximations. q(R) 
is given by the radius of the particle and the method of charging, as can 
be seen from Figure 1. If, however, the spread of particle size is 
narrow enough, then an average q per particle may be used as follows: 



I (go) 

g 

v(«u) 



= * J 



R'(go) 



n(R)dR 



or 



I g H/in p 1 R’(o)), 



(5) 



v(cjj 



= - Jo “(b)® 
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Therefore, the ratio of the "adjusted" value of the current to the 
fluid velocity represents the integral of the particle size distribution. 
When the radius R'(uu), equivalent to a particular frequency in question, 
is larger than that of the particles, the ratio of Equation 5 will be a 
constant. When it lies within the range of the particle size distribt- 
tion, however, a decreasing value for this ratio will be seen as the 
frequency increases. At a high enough frequency, corresponding to radii 
smaller than the particles, the ratio will go to zero. 

The hardware and the measurements are treated next before Equation 5 
is further discussed. 

IV. Experimental Set-Up 

In the EGD apparatus, water droplets are introduced into an air- 

streams as depicted in Figure 3* The injector consists of an aerosol 

flow in conjunction with an ionizer. Saturated steam is injected through 

12 

a corona discharge into the turbulent wake of a cylinder. Electrons 
from the needle drift toward the attractor and attach themselves 
to the then-growing water droplets forming negative charged particles. 
This arrangement with the steam as the source of particles permits some 
variation of particle size distribution by mere manipulation of the steam 
pressure and temperature. 

The anemometer is tested in a plexiglass wind-tunnel. The inside 
dimensions are five cm by ten cm in cross section and 25 cm in length. 

The test section contains the steam injector unit. The aerosol injector 
nozzle consists of a stainless steel ring 2.7 mm in diameter. This 
arrangement is shown in Figure 4. 

The steel ring also serves as the corona ring positive electrode 
(attractor) of the aerosol ionizer. The negative electrode consists of 
a stainless steel needle. Typically at five millimeters downstream from 
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FIGURE 3. EGD GENERATOR SCHEMATIC 




FIGURE 4. AEROSOL INJECTOR 



the cylinder nozzle, a third electrode, in the form of a slender rod 
1.5 mm diameter, collects the current. This separation distance is far 
enough downstream for equilibration between the particles and the flow to 
take place and close enough for substantial currents to be collected. 

The air flow is fed through a temperature regulation stage for the 
control of the relative humidity since it was found that water droplets 
were somewhat affected by the humidity of the air during mixing outside 
of the nozzle. The humidity of the air was held at 32$ because results 
were most reproducible at this value. (Humidity control would not be 
necessary with solid particles.) The water-air mixture flows out of the 
test section into the atmosphere. 

The output from either the hot-wire or the EGD-collector probe is 
fed into a black box which houses the circuitry shown in Figure 5- The 
purpose of this black box is primarily for the protection of instrumenta- 
tion downstream of the probe. The spectral data are obtained from a 
frequency analyzer. It is a General Radio Company Type 1920 Real Time 
Analyzer. It consists of a multifilter band-pass section and a multi- 
channel RMS Detector. The signal from the preamplifier is fed into the 
multifilter, separated according to frequency, then fed into the RMS 
detector where the data are converted to binary digital form and the RMS 
level computed by true linear integration. 

Hot-wire units of the constant -temperature type with frequency 
responses over 50 KHz were utilized. The hot-wire readings were used as 
a standard of comparison for the EGD measurements but, because the hot- 
wire calibration is not appropriate for two-phase flows, air was 
substituted for the steam in the cylinder. The air flow was matched so 
as to give the same momentum as the steam. The absence of droplets also 
helped prolong the life of the hot wires which was very marginal to begin 




FIGURE 5. 

ELECTRICAL SCHEMATIC OF COLLECTOR 



with in the intense vorticity of the wake of the cylinder. The 
resulting hot-wire profile, after the temperatures of the two 
jets were equalized, is shown in Figure 6 for both the main air and 
nozzle air. There are two dominant frequencies apparent on this 
figure, namely, 1.2 KHz and lOKHz. The first one has the value of 
twice the Strouhal frequency for the cylinder and is appropriately 
so because measurements were taken along the centerline of the 
wake where the vortices shed-off by the cylinder on alternating sides 
add up. The Reynolds number based on the cylinder diameter for this 
flow (56 m/sec) is 8 xllO^ so that the boundary layer on the cylinder 
is laminar before separation. The second important frequency in Figure 6 
is associated with the nozzle for the steam- -the inner jet. The 10 KHz 
peak appeared whether or not there was air issuing from this nozzle. The 
exit diameter for this nozzle is 2 mm which is some ten times smaller than 
the cylinder diameter and appropriately the frequency is nearly ten times 
higher. It was observed that 10 KHz is a high enough frequency so that 
the dissipation of the vortices is high; at a station 10 mm away from 
the nozzle exit the 10 KHz peak disappears whereas the 1.2 KHz peak 
remains nearly as intense. 

V. Results and Discussion 

The operation of the EGD- anemometer was investigated for frequency 

response and for current convection efficiency. The high frequency response 

was most noticeable when the steam was very close to the saturated 

condition before expansion (charges of 1°C noticeably affected the results). 

Because of this, it was decided to run the steam at two stagnation 

pressures to achieve different size of particles. Previous work in our 

13 

laboratory indicated that for the aerosol injector used and for an air 
flow of 56 m/sec, smaller droplets appeared at steam reservoir pressures 
below 1.6 atm. 



16 




FIGURE 6. COMPOSITE OF HOT-WIRE ANEMOMETER 

DATA 
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The method of comparison between the EGD and hot-wire data is based 



on Equation 5* 



T 

g 



(<ju) is defined as an RMS parameter, namely. 



g 



I /An q 
(CD) = -gP-— P- 



n 



where 



I = spectral probe current 

A = collector area 
n = total number of charge particles per unit volume 

q = average charge per particle 
v n = flow velocity (RI>G) 

Letting u be a reference voltage (one millivolt in this case), one 

obtains for the value of T in dbs , 

g ’ 




I /An q v 

20 log.. n [T (to )] = 20 log -S2 2_ _ 20 log,. — 

to 10 g v n' & 10 u & 10 u 



— ^ ^ 



J V- 



(6) 



T (to ) 
g v n / 



W 



T. (io ) 
h v n y 



where 

TgC^n) = measurement of the EDG-probe in dbs 

T, (co ) = measurement of the hot-wire in dbs 
h' n y 

The points shown in Figure 7 through 10, have the original EGD-probe 
data and a curve corresponding to the difference between those data and 
the hot-wire values of Figure 6. Twice the Strouhal frequency was chosen 
as a reference point in this comparison because it is the most intense 
point and because frequencies higher than 1.2 KHz are most sensitive 
to changes in the steam parameters. 
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FIGURE 8. COMPARISON OF EGD AND HOT-WIRE 

DATA 
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Figures 7 and 8 show typical results for a steam reservoir pressure 
of 1.85 atm but with differing temperatures and corona voltages. At 
1.85 atm, the corresponding saturation temperature is ll8°C whereas the 
measured values were U6°C and 1T7°C for the data on Figures 7 and 8 
respectively. Figures 9 and 10 show typical results at 1.55 atm. Here 
the saturation temperature is 113°C so it is clear that some super- 
saturation was achieved at both operating pressures before the expansion. 

If the curve for T were a horizontal line, then the two sets of measure- 

g 

ments would only differ by a constant and the EGD anemometer would 
essentially reproduce the hot-wire anemometer measurements. An examina- 
tion of Figures 7 through 10 shows that the present EGD-probe can measure 
intense turbulent fields with some degree of success. Figure 8 shows a 
correspondence of within + 2dbs up to 630 Hz and Figure 9 shows a 
correspondence within + 6dbs up to 20 KHz. The EGD -probe's response can 
be improved by operating with droplets of considerably smaller size 
distribution as will be discussed below. 

The technique outlined earlier in this paper is entirely feasible 

for the measurement of particle size distribution. It is evident from 

Figures 7 and 8 that the data at I .85 atm are well represented by the 

assumptions stated in Section III. There it was said that when the 

particle size distribution is such that it overlaps into the region of 

interest in the frequency scale (from Equation 5) then the amplitude of 

the signal will diminish in a monotonic fashion. The data in Figure 7 

may be interpreted as indicating a distribution with center radius of 
“5 

2 x 10 m and with a half -width of about 1.5 x 10 m (assuming a 
Gaussian distribution for the charged droplets and a corresponding 
average charge per droplet). The data in this figure have a correspondence 
with the hot-wire only up to about 80 Hz, which means that R'(u)) is 
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2.9 x 10 ^ m. R'(uj) is essentially the higher value of the radius at 
the half -width. The data from Figure 8 has a correspondence with the 
hot-wire up to 630 Hz and this yields and R'(a>) of 1.2 x 10 - ^ m. The 
data at 1.55 atm, while having a better overall frequency response, do 
not show a monotonic decrease with increasing frequency. This implies 
the presence Q f more than one Gaussian distribution for the particles 
(the possibility of droplet fragmentation in the wake of the cylinder has 
been proposed '"’). It is evident that conditions outside of the nozzle 
strongly affect the droplet size and the statistical sample of data 
available at the present is too sparse to allow further conclusions. It 
can be seen from Figure 2, however, that a more favorable particle size 
distribution (i.e., centered about 10 m) would enhance the EGD- probe’s 
response to match that of the hot-wire up to 20 KHz. 

Turning now to other operating characteristics of the anemometer, 

Figure 11 shows the effect of the density of the steam on the current 

convection efficiency (I / I ) for the two-phase fluid. In addition to 

§ * 

the sensitivity of this efficiency to small changes in voltage, the effect 

of an increase in density is to require a higher voltage for the same 

efficiency. This shifting of the curves to the right is similar to that 

14 

reported by Whitby for the effect of pressure (temperature was held 
constant) in a single-phase fluid (air). Steam saturation effects are 
not apparent when pressure and temperature data are converted to density 
but they play an important role . 

Figure 12 shows that a definite optimum voltage exists for operation 
with a good high frequency response. After a "triggering" voltage is 
attained the collector current increases linearly until finally a plateau 
is reached. The graph shows that the optimum current efficiency does not 
occur at optimum voltage although some preliminary observations indicated 
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FIGURE 12. OPTIMUM VOLTAGE AND CURRENT 



that acceptable results occurred at a current convection efficiency of 
0.2 to 0,5 with higher collector currents. Figure 13 shows the effect 
of efficiency on the high frequency response for three different 
temperatures with other parameters held constant. 
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CURRENT EFFICIENCY 



Appendix A 



One -Dimensional, Spectral Analysis 

For the case where the energy exchange between the flow and the 
electric field is slight, only the particle momentum equation is 
considered. The one -dimensional case depicted in Figure l4 is a useful 
simplification--it represents an extreme case where the entire fluid 
velocity is opposite to the electric field. In turbulent flows, this 
would be true only for part of the time. But one is interested here on 
an estimate of the maximum slip and hence on a description of the 
maximum interaction. The particle analyzed is assumed to be 
large enough to experience fluid dynamic drag forces (i.e., its diameter 
is greater than the mean free path of the gas molecules) but small 
enough so that pressure forces may be neglected. Stoke 's drag is assumed 
applicable even though this overestimates the actual drag for Reynolds 
numbers (based on the particle diameter) which exceed one. In a restricted 
sense the particles size will be allowed to range to below the mean 
free path of the gas'*’^. 

The free stream velocity "v" is divided in to spectral components at 
preselected frequencies u>^, ... o> n as follows: Let the velocity 

at a given location be 



v(t) = EV(o> )f(u> ,t) 

where 

V(cu ) denotes the amplitude 
and f(u>^,t) is a periodic function with 



|f(«> n ,t)| = 1 '° 
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Constant area channel 




FIGURE 14. ONE-DIMENSIONAL MODEL 
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Now, at a given frequency the particle momentum equation can be 
written as 

dv (u) ,t) 

m p ~ " ' P ' dt" = 6nRpJv(qj n ,t ) ^ - v p (w n ,t)]^ - ^qE^ (Al) 

Stoke 's Drag Electric Force 

nip = mass of the particle 

R = radius of the particle (spherical shape) 
p ro = viscosity of the fluid 

The electrical force is assumed steady here. There may he space charge 
fluctuations which influence "E" locally but these are assumed to be 
outside of the frequency range of interest. 

The particle velocity can be similarly written as 



v = E V (to ) g (to ,t) 
p n p v n /BV n’ ' 



where V is the amplitude and again g(to n ,t) 
ls( ci) n ,t)| = 1.0 so that Equation Al becomes 



is periodic with 



where 



m V (to ) 
P P 11 



dg(co t) m 

■ = — [v(to )f(co ,t) - V (to )g(co ,t)] 
r v n' v n’ ' v n /e>v n’ ' 
v 



dt 



qE 

(A2) 



T 

V 




(A3) 



which is the equilibration time for a constant velocity and no electric 
fields . 

Now, if the particle mass (m ) is vanishingly small, one gets the 

P 

usual definition for the mobility from Equation Al, namely 



D 



= v - 



V 

p 



= pE 



(A4) 
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where 



Vp = drift velocity 



(j, = particle mobility = g - ^- 



Before considering the inertia term as a perturbation, Equation A2 
will be nondimens ionali zed. Let 

A 

t = tco 
n n 

so that for the u> particle, in simplified notation, one gets 



dg m 

m V u> -rAr = r^-[v - v ]-qE 

p p n dt t n pn 



n v 



where 



or 



v = v(u) ,t) and v = v (cu ,t) 
n ' n ’ pn p v n ’ ' 



m V <a> v dg m 

(. P. P n\^n _^ [v . v ] _ ! 

\ qE /dt t qE L n pn J 



qE /dt t qE 
u n v* 

0(e) 0(1) 



(A6) 



0 ( 1 ) 



From the above equation, one can see the perturbation parameter 
explicitly, namely, 

m V co 

c - -BJ 2-E 

- qE 

Now let 

- e„ (0) + £ ®n (l> + °< c2 ) 

Substituting into Equation A6, collecting terms, and neglecting those of 
2 

0(e ) 5 one obtains 



m 



[v - V g (o) ] -1=0 



t v < 3.E n pn n 



(A7) 
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and 



ar [e n <0)] - - iVv. (l) J <*> 

n v 

From the second of the two equations (Equation A8) 

(1) V E d ®n <0) 

v - - m”v ~a$— 

p pn n 

Because E is assumed steady, Equation A7 replaces Equation A5, 
which defined the mobility. Moreover, 



V 

pn 



dg 



(o) 



n 

d$" 



dv df 

_ 22. = V -^2. 
dt n dt 

n n 



so that Equation A8 becomes 



and 



g. 



( 1 ) = _ 



t qE V df 
v 1 n n 



m „2 
P V 
pn 



dt 



n 



pn 



V g = V Cg + eg + 0(e 2 )] 
pn°n pn n v 7 



or 



v 

n 



T q E 
m 



V df 
n n 

V etT 
P n 



m V (jo _ 

p p n . 
qE m 
P 



pn 



which is equivalent to 



= v 



df 

M-E - T OJ v 

v n n dt 



n 



n 




The above equation is consistent with steady flow results and represents 
a first order perturbation of the particle velocity due to inertia effects. 
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Now, from Equation A4 the spectral drift velocity becomes 



Vv*) ' 



V f - 
n n 



V g 
pn°n 



df 

= [iE + t a) V 

v n n at 

n 



The spectral slip then becomes 



pE 
'n ~ V 



df 



n 



S = + V, ar 

n 



But by construction, it is reasonable to assume that 
df 

lari - 1-0 



hence 

pK 

S I = + T U) 

n'max V v n 
n 

where E = is the breakdown field strength. Now is the amplitude 

at the given frequency of the main flow and it is understood that where 

V n is not significant (i.e., zero) the slip parameter has no meaning. 

Now p and t are a function of the radius and p in addition of the 
v 

charge per particle which in turn is a function of the radius and the 
method of charging. 
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Appendix B 

Effect of Particle Size Distribution 
The current arriving at the probe is given by Equation 3> assuming 
one -dimensional flow. 



I g H = Jo v p KR)q(R)n(R)dR 



(3) 



where 



J n(B)aR ■ n 



and v p is the RMS value of the particle velocity. 
In Appendix A, it was found that 

dv 



v = v - llE - t — — 
pn n ^ v dt 

Since |J.E is independent of time, under the present assumptions the 
unsteady part of the above equation is 

I 



dv 



V = Y 

pn n 



n 



v dt 

'where the prime denotes the fluctuating components of the velocity. 
Therefore > 



(A9) 



-2* « 1 - r i(inv') 
v v dt ' n 

n 



(Bl) 



l6 

Now from a linearized version of the Navier-Stokes Equation , 



v = V exp(-icu t) 
n n ' n 7 



where 



i = imaginary index 



Thus , 



-^-(An v ) = -ioj 
dt ' n n 
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and Equation B1 becomes 



= i + iu) t 



n 



n v 



Note that the above is time independent, 
taking the real part only, is 



is = (1 - o> 2 T 2 )i 

v v n v 



so that its RMS value , 



(B2) 



n 

Figure 15 shows a plot of this function for two arbitrary frequencies. 

It can be seen that because of the strong dependence of on the radius 
the function plots as a low-pass filter which may be idealized for the 
present purposes to a step function. 
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FIGURE 15 . PLOT OF (l-rjtu 2 ) 1 ' 2 VERSUS RADIUS FOR 500 HZ AND 16000 Hz 
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